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The spectral properties of a tractable collective model Hamiltonian are studied. The potential 
energy is truncated up to quartic terms in the quadrupole deformation variables, incorporating 
vibrational, 7-independent rotational and axially deformed rotational structures. These physically 
significant limits are analysed in detail and confronted with well-established approximation schemes. 
Furthermore, transitional Hamiltonians in between the limits are presented and discussed. All results 
are obtained within a recently presented Cartan-Weyl based framework to calculate 5(7(1, 1) x SO{5) 
embedded quadrupole collective observables. 
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I. INTRODUCTION 

The low-lying nuclear structure properties of heavy 
and medium-heavy atomic nuclei away from the shell 
closures are dominated by collective modes of motion 
P, ^ . Although contemporary nuclear shell-model calcu- 
lations have pushed the barriers of computational feasi- 
bility with the advent of large-scale shell-model schemes 
Q , this region of the nuclear chart is still unaccessed ter- 
ritory for conventional shell-model calculations. This is 
mainly due to the large degree of collectivity, requiring 
an exceeding amount of shell-model configurations with 
respect to the convergence of the ground- and first ex- 
cited states. To describe the structure of these atomic 
nuclei, one can use (relativistic) mean-field based models 
[1, Si 0, Si or truncated shell-model descriptions, such 
as e.g. the symmetry-based Interacting Boson Model 
(IBM)[|. These models are remarkably successful in de- 
scribing a diversity of collective structures, such as vibra- 
tional and rotational bands or beyond (triaxiality, shape 
coexistence, ...). 

The variety of collective structures arising in the IBM 
triggered numerous studies on the geometrical properties 
of its Hamiltonian as a function of the model parameters. 
For this purpose, the coherent state mean-field approx- 
imation of the many-body problem [HI, [13 is very 
well suited as it couples the fundamental s and d bosons 
of the IBM to a quadrupole deformation field a. The 
global minimum of the energy surface constructed from 
this coherent state gives a good approximation of the 
ground-state energy of the Hamiltonian and, in addition, 
we obtain a geometrical interpretation in terms of the 
quadrupole variables a at the minimum. 

Eventually, the semi-classical limit approach of the co- 
herent state formalism ^T§\ fully establishes the connec- 
tion between the mean-field total-energy surfaces of the 
IBM and the potential-energy surfaces of the geomet- 
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rical Bohr-Mottelson model (BMM)[Ti, [H [Tl|. As a 
major result, it has been demonstrated in a series of 
seminal papers [igj that the critical points along 

the symmetry-transition lines of the IBM can approx- 
imately be described by introducing flat potentials in 
the BMM which are analytically solvable, the so-called 
E{5), X{5) and Y{5) solutions. This observation has 
instigated a number of recent theoretical as well as ex- 
perimental studies on the topic of quantum shape phase 
transitions. From a theoretical point of view, large effort 
has been put in the identification of analytically solv- 
able potentials in the BMM [13, [2l| , the investigation of 
the critical-point structure of the total-energy surfaces 
[H m m, Hi, m, [13, mi and the exploration of spectro- 
sco pic prop erties of the IBM along the transition paths 
[H, [11 [30[, [m, [33, [H [H, [H. From the experimental 
side, a vast number of studies have been devoted to the 
signatures of critical behaviour in atomic nuclei and we 
like to refer the reader to Ref. 36] for a timely review. 

Despite being based on different underlying physical 
principles, the IBM and BMM have a very similar al- 
gebraic structure [sj- It is well-established that the 3 
branching limits of the IBM can be associated with lim- 
iting structures in the BMM. The U{5) limit describes 
harmonic quadrupole vibrations around a spherical min- 
imum, the 0(6) limit can be related to 7-independent 
/3-deformed rotational motion and the SU{3) limit has 
a close connection with the rotation-vibration model 
(RVM) limit m, [H of the BMM. 

This paper aims to explore the spectral properties of 
a transparent though inclusive collective Hamiltonian in 
the BMM, capable of reproducing the 3 forementioned 
limits of collective structure. A good candidate for this 
purpose is the tractable collective Hamiltonian [i^, ex- 
pressed as a polynomial up to (a • a)^ in the potential 
energy. The use of the class of polynomial potentials is 
legitimate from a physics point of view since every realis- 
tic potential or interaction can be expressed as a Taylor 
expansion up to a suSicient degree within the region of 
particular interest. The problem then reduces to the con- 
struction of a matrix representation of the Hamiltonian, 
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which can be done makinguse of a harmonic oscillator 
algebra in any dimension |4l|] (see also Sect. [IT]). More- 
over, the approach enables a clear-cut connection with 
the coherent state total-energy surfaces of the IBM. 

The paper is organised as follows. In Sect. |TT1 we reca- 
pitulate the underlying concepts of the collective model 
and its implementation. For this purpose, we will use 
a recently presented method which treats the collective 
model by means of a pure algebraic technique within the 
Cartan-Weyl basis [13, Sll. In Sect, [m] we will dis- 
cuss the tractable collective Hamiltonian with its spectral 
properties for the physically relevant model parameter 
space and present our conclusions in Sect. IIVI 

II. THE COLLECTIVE MODEL 

The collective model is a macroscopic model in the 
sense that it considers the atomic nucleus as a macro- 
scopic object with a well-defined surface, much alike a 
charged liquid drop. It is assumed that the strong attrac- 
tive nucleon-nucleon interaction roughly favours a spher- 
ical shape, so the radius R{9, (j)) of the nucleus can be 
expressed as a multipole expansion around a spherical 
shape 

i?(0,</.)-i?o |^i + ^«i-r'(e,0)^ , (1) 

with i?o the mean radius, Yl^{9^ (f) the spherical harmon- 
ics and aim the collective variables of multipole order 
I and projection m. The dot denotes invariant angu- 
lar momentum coupling ai ■ bi — {—YV^l + l[ai&;]Q°''. 
Truncating the multipole expansion for the radius ([T]) 
at quadrupole order, we obtain the equation of an el- 
lipsoid, as long as the quadrupole variables a2/j are suf- 
ficiently small (we will abbreviate a2/i by from this 
point). The monopole variable aoo is fixed by imposing 
volume conservation and the dipole variables ai^ merely 
describe small shifts of the centre of mass. From a quan- 
tum mechanical point of view, eq. ([1]) is to be interpreted 
as a dynamic rather than a static deformation of the 
atomic nucleus, the dynamics being determined by the 
Bohr Hamiltonian [l^, [l^ 

HB = f + V{a), (2) 

The quadrupole deformation around the spherical shape 
is assumed to be small, justifying the S'0(3)-scalar Taylor 
expansion of the potential around = 

V{a) — C2{a ■ a) + C3,{[aa\'^ ■ a) + C4(a • a)^ 

-I- C5([aQ;]^ ■ a){a ■ a) + CQ{a ■ of' -f . . . , (3) 

as has been proposed in the General Collective Model 
(GCM) by the Frankfurt group [S H E S IH , which 
can be regarded as an extension of the BMM. This exten- 
sion is most manifest in the description of the kinetical 
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f - 5i^(7r • tt) + Bz{[%af ■ tt + h.c.) + . . . , (4) 

with the canonic momenta tt defined by means of the 
standard Heisenberg-Weyl commutation relations (we 
will omit the operator sign not to overload the notation) 

[tt^, a^] = -iMp^, [tt^, 7r^]=0, [a^, ai.] = 0. (5) 

In the present work, we concentrate on a truncated ver- 
sion of the collective Hamiltonian, as has been discussed 
by Caprio f4(|. Here, the potential energy ^ consists 
of terms up to (a • a)^, which is sufficient to describe vi- 
brational, 7-independent and axially deformed rotational 
structures. The kinetic energy term ^ is chosen in line 
with the standard quadratic expression, as proposed by 
Bohr and Mottelson 

In general, the Hamiltonian ([2]) is not analytically solv- 
able, so a suitable basis is required for numerical diago- 
nalisation. The five dimensional (5D) harmonic oscilla- 
tor basis functions provide a genuine set of basis func- 
tions for this purpose as they can be generated by means 
of an su(l, 1) X so(5) algebra, which is tailor-made for 
polynomial Hamiltonian problems in an N dimensional 
Euclidean space [4l|. However, due to the dimensional- 
ity of the problem, the construction of an orthonormal 
set of 5D harmonic oscillator wavefunctions is a non- 
trivial and even ambiguous task if the quantum number 
L of the angular momentum algebra so(3) is to be pre- 
served. Therefore, a number of different techniques have 
been proposed and extensively discussed in the literature. 
The basis functions up to L = 6 have been determined 
first by Bes by means of a coupled differential equation 
technique (48j . Later on, the importance of the underly- 
ing algebraic structure of the 5D harmonic oscillator has 
been fully appreciated, leading towards the construction 
of explicit basis wavefunctions from basic tensorial build- 
ing blocks [i^, [13, [in or from a coherent state formalism 
[52 | , with good orthonormality properties ^53', '53| . Also 
the vector coherent state formalism 55, 56, 57, 58| and, 
more recently, the algebraic tractable model 0, [60, [Ml 
are suited for the construction of quadrupole harmonic 
oscillator representations with good angular momentum. 
Alternatively, we have shown recently (4^ . [4^ that it is 
possible to calculate the necessary matrix elements of 
all the physical observables with a pure algebraic tech- 
nique i.e. without the need for explicit basis wavefunc- 
tions. This technique transforms the angular momentum 
coupled su(l,l) X so(5) algebra (see eqs. [6]|9l) into the 
Cartan-Weyl basis [62, [SJ, [6J| to calculate the matrix el- 
ements by means of an intermediate-state method. The 
si^ l, 1^) generators are given by the standard expressions 

Bt^ifot.fot, B = fb-b, Bo - i(6t -fo + b-fot), (6) 

with 6t and 6^ the bosonic creation and annihilation op- 
erators 
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and k the spring constant of the 5D harmonic oscillator. 
The Cartan-Weyl realisation of the so(5) algebra is given 
by 112,141 



X 



KV2L+1+V30+1), Y+ 



-'^0 = To{Lo + 30o), Yq = jq{3Lq — Oq), 



4^0+3, 
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with the and respectively the angular momentum 
and octupole generators 



= VTO[b^b] 
= VTO[b^b] 



(1) ^ 



(3) 
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From this definition it can be seen that the sets of 
generators {Xo,X±} and {lo,i^±} span the so(4) = 
su{2)x X su(2)y algebra, which leads to the following 
group reduction 

SU{1, 1) X SO{5) D 50(4) ^ SU{2)x x SU{2)y ■ (10) 

(n,v) V X (X,Mx) (XMy) 

The Cartan-Weyl basis \nvXMxMY) of the s'u(l,l) x 
so(5) algebra is then defined by [isj 



BolnvXMxMy) = i(2n + t; - 
C2M5)]|nuXMjfMy) = w(u + 3)|nwXMxMY), 
C2[sM(2)x,y]|nuXAfjfMy) = X{X + l)\nvXMxMY), 
XolnvXAIxMy) = MxlnvXMxMy), 
YolnvXAIxMy) = MyWvXMxMy) , (11) 

with {n,v) e fe and X = 0,i,l,...,f. {Mx,My} 
follow the standard SU{2) reduction rules with respect 
to X. It is convenient to calculate the matrix elements 
of and tt^ (or equivalent 6^ and 5^) in this basis, as 
both operators carry good bitensorial properties within 
the SU{2)x X SU{2)y reduction according to Racah's 
definition ^67i] . From this point, we will proceed with a^, 
though the results are generally valid for tt*, bj^ and b^i. 



[X±,a^^t] - V(ATM)(A±M+l)a^±i., 



[l±,a^^] = ^/(AT^^)(A±z. + l)ar:^±l, 



(12) 
(13) 
(14) 
(15) 



where the collective variables have been relabelled as 
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which points out that the 5 projections of can be sub- 
divided into the 4 components of a {-kh} bispinor and a 



single {00} biscalar. These bitensorial properties facil- 
itate the calculation of the matrix elements of the col- 
lective variables (and canonic conjugate momenta) con- 
siderably since we can revert to double-reduced matrix 
elements making use twice of the Wigner-Eckart theo- 
rem 



{nvXMxMY\a^t\n'vX'MxMY) 
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with (I) ^ 2X ~ Mx - My and {nvX\\\a^\\\n'v' X') the 
double-reduced matrix element. Explicit analytic expres- 
sions for these double-reduced matrix elements can be ob- 
tained with an intermediate-state method and for further 
details of this derivations, we refer the reader to [12, E^- 

The main difference between previously proposed 
methods and the present Cartan-Weyl based method is 
the embedding of the angular momentum algebra so(3). 
Whereas in the previous methods, explicit wavefunctions 
carrying good angular momentum are constructed to cal- 
culate the matrix elements, this step has been withdrawn 
in the Cartan-Weyl scheme. So, regardless the unphysi- 
cal nature of the quantum numbers X ,Mx and My ^ the 
weight basis is unambiguously defined and leads towards 
an algebraically and computationally enhanced calcula- 
tion of the matrix elements (see eq. (HZ])). Once the ma- 
trix elements are computed, the Cartan-Weyl basis needs 
to be transformed back to the angular momentum basis, 
which can be done by diagonalising the operator 



L • i = 4^2 - 3[(Xo - 3^0 + \){X^ 
-h4V3[r_i iX_ + Ti_iX+\ 



l2Ti 



lT_i 1. 

I 2 2 



(18) 



After the rotation to the physical basis, the diagonali- 
sation of the Hamiltonian proceeds within the separate 
angular momentum L subspaces, similar to previously 
proposed methods. 

The non-compactness of the SU{\, 1) underlying sym- 
metry of the BMM requires extra care with respect to 
the diagonalisation of the Hamiltonian. In principle, the 
Hamiltonian lives within an infinite dimensional Hilbert 
space, so we need to ascertain that the eigenvalues in the 
restricted space are sufficiently close to the exact values in 
the infinite space. Therefore, we gradually enlarged the 
Hilbert space with subsequent harmonic oscillator shells 
until convergency is reached. To find an optimised value 
for the basis harmonic oscillator spring constant, we ap- 
plied the method of Margetan and Williams [1^ . For all 
the calculations in this work, this resulted in a conver- 
gency of 1 eV for all depicted states within 100 harmonic 
oscillator shells or less. 
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C2 



y^phor 200 MeV MeV MeV 
Ky.ind -200 MeV MeV 2500 MeV 
Kot -200 MeV 700 MeV 2500 MeV 



TABLE I: Parameters used in the potential (|20p for the 3 

hmiting C8.S6S V^pher; V^_ind 3^11(1 V^ot • 



III. THE MODEL SPACE OF THE 
COLLECTIVE MODEL 



We intend to study the spectral properties of the col- 
lective Hamiltonian 



(19) 



Expressed in the intrinsic framework, the potential en- 
ergy Via) is written as 



(20) 



V{P, 7) = - V fCs^ cos 37 + C4/3 



which clearly points out that this potential can cover vi- 
brational, 7-independent and axial rotational structures 
for different values of the parameters {c2, C3, C4}. We will 
refer to these potentials respectively as the V^phcr, V^7-ind 
and Vrot hmiting potentials of the collective model. How- 
ever, it should be emphasised that these limits are not to 
be considered as genuine branching limits associated with 
a symmetry algebra such as is the case e.g. in the IBM. 
The algebra su(l, 1) x so(5) in this case is an underlying, 
rather than a spectrum generating algebra, except in the 
case of the 50(5) invariant 7-independent potentials of 
the collective model (eq. (I^Ul) with C3 = 0). The specific 
choices for the parameters of the limiting potentials are 
presented in Table HI The value of f? IB2 is chosen as 4 
keV for all calculations, which agrees with values already 
used in realistic calculations [69|. The motivation for 
these parameter-choices is twofold. Within the restric- 
tions of the given limit under study, the parameters have 
been chosen equal wherever possible. This allows us to re- 
late the differences arising in the structure to the particu- 
lar parameter which has been varied. On the other hand, 
although the sets of parameters give rise to schematic 
potentials, the connection with experimental observables 
was never neglected. As a consequence, the three Hamil- 
tonians can each act as a starting point for a profound 
study of atomic nuclei where typical fingerprints of the 
given limiting cases have been observed. Though the pa- 
rameters {c2, 03,04} might appear to be rather large at 
first sight, one needs to realise that the collective coor- 
dinate /? describes small deformations. Therefore, the 
parameters need to be significantly large to contribute to 
the general structure of the potential. 



A. The three limiting cases 

1. Harmonic quadrupole oscillator 

The harmonic oscillator is conceptually and computa- 
tionally the simplest limit. The Hamiltonian 



sphcr 



2B2 



Vs 



sphcr ; 



(21) 



reduces to the generator Bq of the su{l, 1) algebra, which 
results in the well-known linear spectrum, depicted in 
Fig. [H Apart from the spectrum, all B{E2;Li Lf) 
transition values (relative to B{E2;2i Oi)) are pre- 
sented up to harmonic oscillator shell n = 3. For n — 4, 
only those states that can be organised into a band are 
shown (where a band is defined by following the cascade 
of large relative B{E2) values). For this reason, only 
the largest B{E2) value is given for this set of n = 4 
states. The excitation energy of the 2i state is 1264.91 
keV and the B{E2]2i Oi) = 0.00158 in units relative 
to (3^ei?g)^/(47r)'^. As all mentioned states through- 
out this paper have positive parity, the parity sign will 
not be explicitly denoted. The energy spectrum can be 
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FIG. 1: The level scheme and B{E2) values of a harmonic 
oscillator potential (I21|l . Energy eigenvalues are given rela- 
tive to the first L = 2 state and B(E2) values relative to 
B{E2-2i Oi). 

determined solely using the su(l, 1) part of the algebra, 
whereas the calculation of the B{E2) values also requires 
the so(5) part. This can be seen from the definition of 
the B{E2) reduced transition probability 

B{E2;L, -^ Lf) = e^ J2 \{LfMf\fiE2)^\UAU)f, 

(22) 

with 



f{E2)^ = ^a;. 



(23) 
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as Qf^ is a v = 1 S'0(5)-tensor connecting different 5*0(5) 
representations with Av = 1. Z is the proton number 
and e denotes the elementary charge of the proton. 

Therefore, as the 5D harmonic osciUator has aheady 
been discussed numerously in the literature, it provides 
a reliable test for the matrix elements in this Cartan- 
Weyl reduction scheme. As a result, we can proceed with 
confidence studying structures which can no longer be 
expressed as a genuine symmetry limit and for which a 
numerical treatment is required. 



2. 'y -independent rotor 

The second limit to be considered is the 7-independent 
rotor, described by the potential V^-y-ind 



23. 



(24) 



The occurrence of a global minimum at non-zero defor- 
mation /3o = makes this potential essentially different 
from the harmonic oscillator potential. Consequently, 
this potential is able to generate solutions which can be 
associated with definite deformations. 

Although the Hamiltonian does not correspond to an 
algebraic solvable limit, still it exhibits some remarkable 
symmetry properties. Similar to the case of the harmonic 
oscillator, the Hamiltonian is an SO{b) invariant, lead- 
ing towards degeneracies within a given representation v. 
This is illustrated in Fig. [2] where for all bands except the 
one in the middle, the states with equal seniority v have 
the same energy eigenvalue. The main difference between 
H^.ind and i/spher with respect to the spectrum generat- 
ing algebra sm(1, 1) is that, although the Hamiltonian can 
entirely be written as a function of the generators of this 
group, it does no longer constitute a dynamical symme- 
try for the Hamiltonian. In this particular case, it causes 
the n ^ states of the harmonic oscillator to be lifted in 
the spectrum. This is illustrated in Fig. [5] by the band 
in the middle, which lies much higher in the excitation 
spectrum than the corresponding band of the harmonic 
oscillator (Fig. [T]). 

Similar to Fig.[l] all eigenstates with an excitation en- 
ergy lower than or equal to the excitation energy of the 
61 state are depicted in Fig. [21 together with all non- 
vanishing B{E2) values. For the higher- lying excited 
states, only those states are given that fit into a given 
band, as well as the corresponding intraband B{E2) 
values to justify this classification. All eigenvalues are 
again given relative to the 2i state, which has an abso- 
lute excitation energy of 233.8 keV, and the B{E2; 2\ 
0i)=0.00759, measured in units i^iZ eRlf j {'^-nf . 

To validate the numerical calculations and to inter- 
prete the results, it is instructive to compare them with 
analytically solvable approximation schemes. Such a 
scheme is provided by the displaced harmonic oscillator 
of the Wilets & Jean (WJ) class of /? — 7 decoupled po- 
tentials [zOl- We can perform a Taylor expansion of the 
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FIG. 2: The level scheme and B{E2) values of a 7- 
independent rotor potential (|24|) . Energy eigenvalues are 
given relative to the first L — 2 state and B{E2) values rela- 
tive to B{E2;2i Oi). 



potential V^-y-ind around the minimum /3o = 0.2 of the 
potential, leading to a WJ displaced harmonic oscillator 
as long as the expansion is truncated up to second order. 
The eigenvalues of this displaced harmonic oscillator are 
given by 



E„ 



hn[u;'{n0 + l) + l{(]'„iv)-Po)i2f3'oiv)-Po)], (25) 



w ith n = ^-4c2/B2, w' = V4-3/3o//?o(«), k = 
^/—4:C2B2/fi, and /3q(w) the minimum of an effective po- 
tential, determined by the solution of 



(26) 



Inserting the specific parameters at hand, we obtain an 
excitation energy of 210.3 keV for the 2i state ((ri/j, v) = 
(0, 1)), and the ratio £;ox(4i)/£^ex(2i) = 2.45, which com- 
pares well with the numerical results (see Fig. [2]). More 
interesting are the predictions for the first excited u = 
state {{npjv) = (1,0)). The excitation energy of this 
state is predicted as 1805.1 keV, which is in reasonable 
agreement with the value of 1700.6 keV of the O3 state 
depicted in Fig. O Here, we need to take into account 
that at this energy scale, we are at the limits of the ap- 
proximation's validity with -4 MeV being the minimum 
of the original potential V-y-ind- Similar conclusions can 
be drawn for the other states in the band built upon the 
O3 state in Fig. [2] and therefore, we refer to this band 
as the /3-vibrational band of the 7-independent model, 
since we can approximately associate the /3-vibrational 
quantum number = 1 to this band. 
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3. Axially deformed rotor 

The third hmit describes axially deformed rotational 
structures. In this particular case, we insert the term 
[aa]^^-* • a in the Hamiltonian 



Hy, 



2B, 



(27) 



breaking all the remaining degeneracies from the 7- 
independent rotor case. Moreover, the classification into 
bands by following cascades of B{E2) values is even more 
pronounced as the bands all occur at different energy 
scales in the spectrum. Whereas the bands built on top of 
the 22 state in the harmonic oscillator and 7-independent 
rotor limit had energy scales comparable with respect 
to the ground-state band energy scale, in this particular 
case, they are observed at much higher excitation en- 
ergies. This is illustrated in Fig. [3] by the dashed box, 
pointing out that the bands on top of the 22 and the O2 
states are higher excited than the 12i state of the ground 
band. 
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FIG. 3: The energy spectrum and B{E2) values of an axial 
rotational potential (|27|) . Energy eigenvalues are given rel- 
ative to the first L = 2 state and B{E2) values relative to 
B{E2;2i ^Oi). 

From an experimental point of view, this might be an 
unsatisfactory situation, as it can occur that the bands 
on top of the 22 and O2 states can appear at a much 
lower excitation energy and come alongside the L = 6 or 
L = 8 members of the ground band. However, the chosen 



potential Vi-ot does not contain many free parameters to 
tune the excitation energy of the associated bands while 
keeping the minimum at a physical deformation (3q ~ 0.1. 
We can as well adjust the mass parameter B2 of the ki- 
netic energy. A decrease of the mass parameter would 
lead to an overall lowering of the excitation energy of 
the excited bands, but this would cause the ground state 
band to lose its rotational character. It should be men- 
tioned that a modified kinetic energy term was proposed 
in the GCM to cope with this issue, leading towards a 
lowering of the excited bands in the spectrum without 
affecting the rotational structure of the bands. However, 
we will not embark into a study of the effects of the mod- 
ified kinetic energy terms in the present work, as we con- 
centrate on the global features of the potentials, rather 
than the phenomenological description of experimental 
data. 

In Fig. [3l only the three lowest bands are depicted, 
with energies relative to the 71.29 keV 2i state and 
B{E2) values relative to B{E2;2i ^ Oi) = 0.01245 in 
units of (3Zei?§)V(47r)2. 

Similar to the 7-independent rotor, the comparison of 
the numerical results with an analytically solvable ap- 
proxmation can shed light on the general structure of 
the solutions. The RVM [sl] is a good candidate for this 
purpose, as it is based on the physical assumption that 
an axially deformed nucleus can be described by means of 
a harmonic oscillator potential, in both 7 and /3 around 
the minimum in the potential V^ot ■ Consequently, we can 
carry out a Taylor expansion of the potential around this 
minimum and compare the approximate solutions with 
the numerical results. These solutions can be cast in the 
well-known expression js^ 



E 



LKn2nQ 



{^\K\ + 1 + 2n2)E^ + {no + \)Ep 



{L{L + l)- K')^e 



(28) 



with K the angular momentum projection quantum num- 
ber along the intrinsic axis and (^2,710) the vibrational 
quantum numbers associated with respectively the vibra- 
tions in the 7- and /3-direction. i?-y, Ep and e are param- 
eters, determined by the shape and localisation of the 
RVM potential. Rewriting the potential V^_rot around 
the minimum (/3o,7o) = (0.264,0) as a Taylor expan- 
sion in j3 and 7, gives rise to the following parameters 
of the RVM: E-^ = 1885.22 keV, Ep = 2093.96 keV 
and e = 19.15 keV. Substitution of these parameters 
in eq. (j28p gives some remarkable results in comparison 
to the solutions of i/iot (see Fig. First, the spac- 
ing of the different levels within a given band is reason- 
ably well reproduced by the parameter e. Second, the 
position of the first K — 2 band is to be expected at 
E^ = 1888.22 keV according to the RVM, which is in 
good agreement with the band built on top of the 22 
state in Fig. [3l More interesting is the classification of 
the first excited K = Q band on top of the O2 state. 
Within the language of the RVM model, we can either 
associate this band with a /3- vibrational (rig = 1 , ^^2 = 0) 
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^ ^ ^ 

FIG. 4: Energy spectra for the transitional Hamiltonians H12 (a), H23 (b) and H31 (c) as a function of the parameter 



or 7-vibrational (no = 0,7i2 = 1) structure. Substi- 
tuting the corresponding quantum numbers (no, n2) into 
eq. (|28|) gives rise to the foUowing predictions: the excita- 
tion energy of the 0^ bandhead of the 7-vibrational band 
is to be expected at 2E^ = 3770.44 keV, whereas the 
0/3 bandhead of the /3-vibrational band can be found at 
Ep — 2093.96 keV. Comparison of these values with the 
excitation energies from the diagonalisation of i?rot leads 
towards the conclusion that the lowest excited K — 
band is best to be associated with the /?-band of the 
RVM. Therefore, we will use the nomenclature of the 
RVM when referring to the different bands in the rota- 
tional limit -ffiot- 

In conclusion, the ability of the quartic truncated col- 
lective model Hamiltonian to cover vibrational, 7- 
independent and rotational limits has been validated 
through the analogy with transparent approximation 
schemes. This enables a clear-cut physical interpreta- 
tion of the bands for each separate limit. Also, it is clear 
from Figs. [T] to [3] that the spectra of the correspond- 
ing limits can be described respectively by the U{5), the 
0(6) and the SU{3) symmetry limits of the IBM, which 
are the well-know vibrational, 7-independent rotor, and 
axial rotor limits of the IBM [7l|. This is not too sur- 
prising as the close relation between the IBM and the 
BMM has long been established (see [131 and references 
therein), mainly based on the general characteristics of 
quadrupole collective models. Intuitively, the connection 
can be made on the level of the potential energy pp]) 



of the BMM. This expression ([^0]) can also be obtained 
as a total-energy surface of the IBM within the coher- 
ent state mean-field formalism 13]. From the point of 
view of the IBM, the ground state properties of a par- 
ticular IBM Hamiltonian are encoded in the minimum 
of these total-energy surfaces, whereas within the BMM 
formalism, the minimum of the potential energy surface 
rather gives the mean value around which the dynamic 
excitations are located. For large values of the mass pa- 
rameter B2 of the BMM, the ground state wavefunction 
is nicely confined within the potential well, with the asso- 
ciated energy eigenvalue approaching the minimum value 
of the potential as the mass parameter reaches infinity. 
In the present calculations, h^/B2 is chosen sufficiently 
small (4 keV) with respect to the typical energy scale of 
the potential (~ MeV). However, we like to stress that 
this is an intuitive argument, rather than a hard proof 
for the manifestation of the connection of the respective 
IBM and BMM limits in Figs. [I]to[31 Presently, we will 
not proceed along this track as the connection is worth 
a study on its own and we refer the reader to recent in- 
depth studies m, [tII on the subject. 



B. Transition paths 

Having determined the structure of the typical limit- 
ing cases in the collective model, it is interesting to see 
how these structures evolve as the potentials gradually 
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ground band 

R{ E2; 2^=> 2|) 



R(E2; 3i=> 2,) 
R( E2; 3,=>2^) 






0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

^ ^ ^ 

FIG. 5: Selected B[E2) values for the transitional Hamiltonians H\2 (a), H23 (b) and H^i (c) along the transition lines, relative 
to the B{E2;2i Oi) value (see equation (l32}). 



change from one limiting case into another. Therefore, 
we construct the following transitional Hamiltonians 



H12 = 




Ki 


^ O^sphor + C^7-ind, 


(29) 


H23 — 




Ki 




(30) 


H31 = 




Ki 


-OK-ot +C^sphor, 


(31) 



for which we calculate the excitation energy and electric 
quadrupole observables along the transition path. The 
results are presented in Figs. IH El [H and [T] For easy 
comparison, the three different transition paths (respec- 
tively H12, H23 and H31) are plotted side by side with the 
same scaling. It should be mentioned that the study of 
H12 is closely related to recent work on phase transitions 
in the Bohr-Mottelson model [t^, where the Hamilto- 
nian is pushed into the domain of large deformations in 
order to clearly identify the critical point in the transi- 
tion and discuss the apparent quasi-dynamical symmetry 
along the transition line. 

We start with the discussion of the energy spectra. In 
Fig. 131 the energy spectra are plotted as a function of 
^. In the left panel (a), the transition from the spheri- 
cal to the 7-independent rotor is presented (see eq. . 
Since 50(5) is a symmetry for both limits, it is also a 
symmetry for every intermediate Hamiltonian H12, which 
makes the seniority quantum number v a good quantum 
number along the whole path. Furthermore, the large 
degeneracy of the harmonic oscillator persists along the 
transition path, except for the band on top of the L — 
state originating from the O2 state of the harmonic oscil- 
lator (see Fig. [1| . One can clearly see from Fig. IHa) how 



this band (depicted in red) decouples from the ground 
band (in blue) to become the /3-vibrational band in the 
7-independent limit. The different nature of this band 
is even more apparent from the pure crossings around 
^ ~ 0.8. Similarly, one can deduce from the relative 
B{E2) values related to this band (see Fig.[5fa)) that this 
/3-vibrational like L = (O2 up to ^ '--^ 0.8) state differs 
structurally from the other one (O3 up to ^ ~ 0.8), as the 
B{E2;02 2i)/B{E2;2i d) value decreases faster 
than the values B{E2;Li ^ {L ~ 2)i)/B{E2;2i d) 
of the ground-band members or even the B{E2; 3i 
"2^2) / B{E2]2i — > Oi). To simplify the notation, it is con- 
venient to introduce the following quantity 



RiE2;L,^Lf) 



B{E2;U^Lf) 



B{E2;2i 



(32) 



Oi)' 

Oi) as plot 



for the B{E2) values relative to B{E2;2i 
ted in Fig. O 

In the middle panel (b) of Fig. [H the transition from 
the 7-independent to axially deformed rotor is plotted 
(see eg. I30p . By introducing the seniority breaking term 
[aa] ■ a in the Hamiltonian, the remaining 50(5) degen- 
eracy of the 7-independent rotor is lifted. This effect 
arises rather promptly with small values of ^, leading to 
the instant development of separate bands in the spec- 
trum. At the right-end of the transition path (^ = 1), a 
clear picture takes on form with distinct rotational and 
vibrational-like bands at different energy scales. The fact 
that the states can now unambiguously be organised into 
bands is also reflected in the R{E2; Li ^ Lf) values in 
the middle panel (b) of Fig. \5\ There it can be seen 
that the value R{E2;22 2i) quickly drops whereas 
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R{E2]3i — > 22) stays reasonably unaffected along the 
transition line, which clearly points out that 22 and 3i 
belong to the same {K = 2) band, depicted in the dashed 
box of Fig. [3l 

In the right panel (c) of Fig. IH we close the circle 
with the transition from the axially deformed limit to 
the vibrational limit. It can be seen from that figure that 
this transition is less gradual, compared to the harmonic 
oscillator to 7-independent rotor transition. In the latter 
case, only the SU (1,1) symmetry was broken, whereas in 
the former case also the SO {5) symmetry is immediately 
broken by adding the [aa]^^^ • a term in the Hamiltonian. 

Before proceeding, there is a peculiarity to be noted 
concerning the identification of the K = Q bands in Fig. [4] 
as either a 7- or /3-vibrational band. From the compar- 
ison of the axially deformed limit with the predictions 
made by RVM, we could unambiguously identify the low- 
est K = (built on the O2 state) as the /3-vibrational 
band. This can further be justified if we follow the tran- 
sition path towards the harmonic oscillator, as depicted 
in the right panel (c) of Fig. U) We notice that, e.g. the 
O2 state of this /3-vibrational band (green lines) evolves 
towards the (n = l,w = 0) SU{1,1) representation of 
the harmonic oscillator, which is basically a /3-vibration. 
However, this is not the only path to reach the harmonic 
oscillator. One could go the opposite way in Fig. [H and 
follow the transition from axially deformed rotor to har- 
monic oscillator limit via the 7-independent rotor limit. 
Following e.g. the O2 state along each transition (middle 
(b) and left panel (a)), we find that in this case, we end 
up in the {n ~ 0,v ~ 3) state, which cannot be asso- 
ciated with /3- vibrations. The solution to this paradox 
lies in the middle panel (b). From this panel, it is clear 
that all states in the 7-independent to axial rotor limit 
are subject to considerable mixing and 'no-crossing' ef- 



fects, due to the large admixture of states with different 
seniority quantum numbers. Moreover, as the O3 state 
of the 7-independent rotor limit clearly evolves from the 
O2 {n = l,v — 0) state of the harmonic oscillator, we 
can state that the O2 and O3 states must have switched 
nature along the transition path from 7-independent to 
axially deformed rotor. In conclusion, the association 
of the K = bands with (3- or 7-vibrational structures 
should be treated with caution as considerable mixing 
effects can arise, perturbing the simple picture of vibra- 
tional motion along the /3- or 7 direction. 

Apart from the relative B{E2) values, it is also inter- 
esting to calculate the absolute B{E2) values, as they 
give an indication of the collective deformation of an 
atomic nucleus [39]. In Fig. [SI the B{E2; 2\ — > Oi) values 
are presented along the transition paths. One can clearly 
see that, going from the harmonic oscillator to the 7- 
rotational limit, the B{E2; 2i ^ Oi) value rises steadily, 
due to the manifestation of a deformed minimum in the 
potential. This rise persists in the transition from the 
7-independent rotor to the rotational limit, as the onset 
of the term [aa] ■ a in the potential breaks the symmetry 
in the 7-direction, driving the minimum in the potential 
towards prolate structures. Finally, along the rotational 
to vibrational transition path, the B{E2; 2\ Oi) drops 
back to the originally value of the spherical harmonic 
oscillator. 

Similar conclusions can be drawn from a study of the 
spectroscopic quadrupole moment of the first excited 2\ 
state, defined by 

g = yi|I(2iA/ = 2|f (£;2)|2iAf = 2), (33) 

where the same linear approximation of T(E2) is used as 
for the calculation of the B{E2) values (see eqs. and 
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([23)) ). The results for Q are presented in Fig.[7l Contrary 
to the other figures, the transition from the harmonic os- 
cillator to the 7-independent rotor is not depicted since, 
within the linear approximation of T{E2), the selection 
rules for a (Aw — ±1) render the quadrupole moments 
identically zero along this whole transition path. Here 
again, we notice that the quadrupole moment sharply 
rises (in absolute value), as soon as the seniority-breaking 
term in the potential is turned on (along both transition 
paths), to reach a maximum at the rotational limit. This 
rise is more pronounced in the transition from the 7- 
independent rotor to rotational limit (b) as compared to 
the transition from the spherical to rotational limit (c). 
This difference stems from the fact that in the former 
transition the potential already exhibits a deformed min- 
imum in the (3 direction, while in the latter transition the 
minimum in the /3 direction is steadily developing from 
the spherical minimum towards a minimum at a distinct 
deformation f3Q. 
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FIG. 7: The quadrup ole mom ents Q of the first excited 2i 
state, given in units ^y^&n/5{3ZeRo)/{i^T). Only the values 
for the transitional Hamiltonians H23 (b) and H31 (c) are pre- 
sented as the quadrupole moments for the harmonic oscillator 
to the 7-independent rotor are identically zero. 



IV. CONCLUSIONS AND OUTLOOK 

In the present paper, we have studied the spectral 
properties of a truncated collective Hamiltonian up to 
quartic terms in the potential energy. This Hamiltonian 
is able to cover three different limits which can be associ- 
ated with vibrational, 7-independent rotational and rota- 
tional structures. This has been demonstrated by com- 
parison of these limits with physically transparent ap- 
proximation schemes and the connection with the three 
branching limits of the IBM. Furthermore, it has been 
studied how the typical fingerprints of collective structure 
(such as e.g. energy spectra, i?(i?2)-values, etc.) evolve 
if the Hamiltonian changes from one particular limit to 
another, which gave rise to some remarkable effects due 
to seniority mixing in the 7-independent rotational to 
axially deformed rotational limit. 

The present results also provide a genuine test for the 
recently developed program which treats the (general) 
collective model within a Cartan-Weyl framework. This 
is necessary if one wants to proceed towards more com- 
plex collective structures, such as e.g. triaxiality and 
shape coexistence, which are prominently present in the 
medium-heavy to heavy isotopes of the nuclear chart. 
Extending the present truncated Hamiltonian to incor- 
porate higher-order terms will be the subject of future 
investigation. 
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